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<^ ' Abstract 

We review recent results on the finite-gap properties of difference operators with 
ehiptic coefficients and give explicit characterization of spectral curves for difference 
analogues of the higher Lame operators. This curve parametrizes double-Bloch so- 
lutions to the difference Lame equation. The curve depends on a positive integer 
number i, related to its genus g hy g = 21, and two continuous parameters: the 
lattice spacing 77 and the modular parameter r. Isospectral deformations of the 
f^ ' difference Lame operator under Volterra flows are also discussed. 

a^ 

o 

o. 

^ 1 Introduction 

^ ■ The spectrum of the Schrodinger operator —d^ + u{x) with a periodic potential u{x) = 

c3 ■ u{x+T) has a band structure: there are stable energy bands separated by gaps. For smooth 

p ! potentials, the width of gaps rapidly decreases as energy becomes higher. However, gaps 

^ I generically occur at arbitrarily high energies, so there are infinitely many of them. 

S^ . Of particular interest are exceptional cases, when for sufficiently high energies there are 

H ! no gaps anymore and their number is therefore finite. Such operators are usually refered 

to as algebraically integrable or finite-gap ones. Their studies go back to classical works of 
the last century. The renewed interest to the theory of finite-gap operators is due to their 
role in constructing quasi-periodic exact solutions to non-linear integrable equations. 

Among examples of the finite-gap operators, the first and the most familiar one is the 
classical Lame operator 

C = ~— + i{i+l)p{x + u'\u,u'), (1.1) 

where p{x\u,uj') is the Weierstrass p-function and £ is a parameter. The potential is a 
doubly-periodic function on the complex plane with periods 2uj and 2a;', where Im (u'/u) = 
r > 0. If a; is real while a;' is pure imaginary, the spectral problem is self-adjoint. The 
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finite gap property of liiglier Lame operators for integer values of I was established in |]T|. 
If £ is a positive integer, then the Lame operator has exactly I gaps in the spectrum. Such 
a remarkable spectral property is a signification of a hidden algebraic symmetry which, in 
its turn, leads to an intimate connection with integrable systems. 

The finite-gap property becomes even more striking for difference operators. A natural 
difference analogue of the Schrodinger equation has the form 

a{x)'^{x + r]) + h{x)'^{x) + c(x)\l'(x - rf) = E'^^x) , (1.2) 

where the parameter rj is the lattice spacing. Let us assume that rj is real, the coefficients 
are real functions of x and c(x) = a{x — if), then the problem is self-adjoint. Let the 
coefficient functions be periodic with a common period T: a{x+T) = a{x), b{x+T) = b{x). 
The difference Schrodinger operators with periodic coefficients exhibit much richer spectral 
properties because the problem has two competing periods (T and 77) rather then one. 
Nevertheless, the class of finite-gap operators survives. 

The structure of the spectrum of a typical difference operator crucially depends on 
whether the ratio T/rj of the two periods is rational or irrational number. In the former 
(commensurate) case one can always set T/r] = Q G Z without loss of generality. Then 
there are no more than Q stable bands in the spectrum. Indeed, set \E'(xo + nr]) = \Ef„, 
a{xo + nr]) = an, etc and rewrite eq. ([L^) in the form 

where an+q = an, bn+g = bn, Cn+q = Cn- Since the coefficients are periodic, one may 
look for solutions in the Bloch form: \E'„ = e^^^^Xn-, where Xn is Q-periodic and k is the 
Bloch momentum. Therefore, the spectral problem is reduced to the eigenvalue problem 
for a hermitian QxQ-matrix. For each real value of the Bloch momentum k the secular 
equation has Q real solutions E = Ei{k). As k sweeps over the Brillouin zone, Ei{k) sweep 
over the stable bands labeled by i. Several neighbouring bands can merge, so the total 
number of stable bands can be less or equal to Q. 

The latter, incommensurate case can be practically realized as a proper limit of the 
former when both numerator and denominator of the fraction T/rj = Q/P tend to infinity. 
The resulting spectra can be (and usually are) extremely complicated chaotic generations 
of Cantor set type. Some of them, like those in the Azbel-Hofstadter problem [^ , though of 
a multifractal nature, nevertheless keep a good deal of hidden regularity revealed in terms 
of string solutions to Bethe equations [^. Very little is known on the spectra of generic 
type; they seem to be completely irregular. In this paper we discuss just the opposite case 
of the utmost regular spectra in the sense that the number of bands is finite although T/t] 
is irrational. Moreover, the number of bands does not really depend on this ratio, being 
determined by another (integer) parameter. The operators with this type of spectra are 
true difference analogues of the finite-gap operators. 

In 0], the following difference analogue of the Lame operator ( |1 . 1|) was proposed: 

i = !l^LfM,* + ^l^±i!l),-*. (1.3) 

Oi{x) 6'i(x) 

Here ^1(0;) = 6i{x\t) is the odd Jacobi ^-function, £ is a non-negative integer. The coeffi- 
cients are periodic functions with period 1. This operator can be made self-adjoint by the 



similarity transformation L -^ g^^{x)Lg{x) with a function g{x) such that g{x + l) = g{x), 
so the spectrum is reaL The operator (pT3|) first appeared in a completely different context 
of representations of the Sklyanin algebra as early as in 1983 0. Namely, L coincides 
with one of the four generators of the Sklyanin algebra in the functional realization found 
by Sklyanin. Remarkably, the operator L for positive integer values of i and arbitrary 
generic r] has 2i + 1 stable bands (and 2i gaps) in the spectrum. The finite-gap property 
of this operator for integer i was proved in |^. It was also shown [^ Q that the Sklyanin 
algebra does provide a natural algebraic framework for analyzing the spectral properties 
of the operator L. (A different algebraic approach to the difference analogues of the Lame 
operators was proposed in 0.) 

Another similarity transformation, L = f^^Lf, where 

fix) = l[9^ix-jv), (1.4) 

i=i 

makes coefficients of the difference operator 

double-periodic functions of x with periods 1 and r. The limit rj —>■ gives the Lame 
operator (p..l|). Indeed, replacing x by x + |r in ( |1.5| ) we obtain: 

L = 2-r]^{C + const) + 0{7]^), r/ ^ , (1.6) 

where the p- function in the C has periods 1 and r (see ([A3|)). 



Let us mention that spectral curves of the classical Lame operator ( |1.1|) and its Treibich- 
Verdier generalizations |^ for small values of £ were studied in p, |TD[. A detailed analysis 
of solutions to the difference Lame equation at £ = 1 was recently carried out in ||ll| . 

The paper is organized as follows. Sec. 2 is a continuation of the Introduction. To draw 
the problem in a broader context, we discuss the general notion of the finite-gap operator. 
In Sec. 3, a family of Bloch eigenf unctions of the operator (|1.5| ) is constructed. These 
eigenfunctions are parametrized by points of the spectral curve. Sec. 4 contains equations 
for the edges of bands and some examples. In Sec. 5 we work out an explicit relation 
between the Bloch multipliers. The form of the result suggests that some hypothetical 
combinatorial identities for "elliptic numbers" may be relevant. At last. Sec. 6 contains 
some remarks on the isospectral deformations of the difference Lame operator. In this case 
the coefficient in ( |1.5| ) has more poles. The location of the poles, however, is not arbitrary: 
they are constrained by locus equations. 



2 A general view of finite-gap operators 

The key idea of the modern approach to spectra of differential or difference operators 
is to regard the solutions \l/(x, -E) to the spectral problem (say, ( |1.2| )) as functions of E 
for any complex values of E and to study their analytic properties in E. In so doing it 
is not necessary to assume that the problem is self-adjoint, so the parameter i] and the 
coefficients may be complex numbers. 



In practice, one may try to construct a family of eigenfunctions \l/ = \l/(x, E,pi,p2, . . .) 
depending, apart from E, on a finite number of additional parameters Pi. For instance, one 
of these could be the Bloch momentum k: \l/(x) = e''''^x{^y E, k), where x{^ + T, E, k) = 
x{x, E, k). Suppose such a family does exist. Then the spectral parameters appear to be 
constrained by some relations Fi{E,pi,p2, . . .) = 0, so that only one of the parameters is 
independent. These relations define a complex curve (a Riemann surface) in the parameter 
space called the spectral curve. The true spectral parameter is a point of the curve. This 
is the proper mathematical formulation of the dispersion law E = E{k). Usually, this 
function is multi-valued. It becomes single-valued on the spectral curve (when the latter 
is well-defined). Moreover, the solution \l/(x,£^) of the spectral problem also becomes a 
single-valued function on the spectral curve. 

In the case of the second-order difference operators the spectral problem ( |1.2| ) has no 
more than two linearly independent solutions. In other words, the function E can take 
any of its values at most twice. The existence of such a function implies that the spectral 
curve is a hyperelliptic curve. Any hyperelliptic curve can be represented in the form a 
two-sheet covering of the complex plane of the variable E: 



.2 



W 



WiE-E,), (2.1) 



where Ei are called branch points. For E = Ei the equation ( |1.2| ) has only one linearly 
independent solution. The curve (|2.1|) is well-defined if the set of branch points is finite. 
Then the curve is algebraic and has finite genus. Equivalently, this means that there exists 
a difference operator W such that W can not be represented as a polynomial function of the 
difference operator in the right hand side of ( |1.2| ), and that commutes with this operator. 
In this case they have a set of common eigenfunctions. Equation (|2.1|) is then lifted to the 
operator relation. The parameter w is the eigenvalue of the operator W on the common 
eigenfunction \l/(x, E): W'^{x,E) = w^!{x,E). 

For self-adjoint spectral problems the branch points Ei are real numbers Ei < E2 < 
E-^ < . . .. In the stable bands the Bloch momentum takes real values. The stable bands 
are segments of the real line [-E2J+I5 -^21+2], i = 0,1, . . ., so the branch points are just edges 
of bands. At the edges of bands, the Bloch solutions \l/(x, Ei) are periodic or anti-periodic. 

The notion of the spectral curve is really useful in the exceptional case of the spectra 
of regular type mentioned in the Introduction. Recall that in the typical case the set of 
branch points may be even uncountable, so the spectral curve, defined as above, does not 
have sense. The very existence of a well-defined spectral curve of finite genus is the precise 
characterization of the finite-gap operators. 

Then it is natural to address the inverse problem: given a hyperelliptic curve of finite 
genus regarded as a spectral curve of some difference operator, to find coefficients of this 
operator, i.e., the functions a{x), b{x) in ( |1.2|) . In this way, one is able to construct a 
representative family of finite-gap operators [112], |T3[. For difference operators this was 



done in |14, O, M. The coefficients of the operators can be expressed through Riemann's 



theta-functions associated with the curve. 

The curve itself does not determine the operator uniquely. There is a remaining finite- 
parametric freedom which can be fixed by some additional data on the curve (essentially, 
a number of marked points). In other words, any finite-gap operator admits a class of 



isospectral deformations. The coefficients of the operator, with respect to the isospectral 
flows, obey certain non-hnear integrable equations. 

Looking for formulas more effective than a bunch of Riemann's theta-fucntions, one 
may inquire whether they can be expressed through simpler functions, for instance elliptic 
ones. For a particular class of curves, namely, for special coverings of elliptic curves, this is 
indeed possible (see e.g. [0). The Riemann theta- function associated with such a curve 
factorizes into a product of Jacobi 6'-functions, so the coefficients of the operator become 
elliptic functions. Moreover, all the family of isospectral deformations of the operator 
enjoys the same property. 

An important example of this phenomenon in the differential set-up is provided by the 
Lame operator ( |1.1|) for £ G Z+ and its isospectral deformations —d1 + u{x) with 

(■((■+1)12 

u{x) = 2 ^ p{x — Xj) + const . (2.2) 

i=i 

The Lame operator itself corresponds to a very degenerate configuration when all the poles 
sit in one and the same point. The isospectral flows are the flows of the KdV hierarchy jT^ 
for the potential u{x). Solving say the KdV equation ii = 6uu' — u"' for u = u{x,t) with 
the initial condition u{x, 0) = i{£+l)p{x — xq), we get a family of Schrodinger operators 
with elliptic potential which have the same spectral curve as the Lame operator. The poles 
Xj (and the constant term) in ( |2.2| ) become t-dependent. By a direct substitution to the 
KdV equation, it can be shown |jl8[ that they are constrained by the conditions 



£{m)/2 
J2 p'{x,-x,) = 0, t = l,2,...,i{i+l)/2, (2.3) 



and obey the differential equations 

i{m)/2 



X 



J 



-12 Y. Pi^i-^k). (2.4) 



k=l,^j 



Eqs. (|27^ ) are the famous equations defining the equilibrium locus of the elliptic Calogero- 
Moser system of particles. From the general theory which connects the pole dynamics 
of elliptic solutions of non-linear integrable equations with systems of Calogero-Moser 
type [|1^, it follows that the connected component of the locus is parametrized by the 
Jacobian of the spectral curve of the Lame operator. So, it is an ^-dimensional submanifold, 
spanned by higher Calogero-Moser flows, in the |£(£-|-l)-dimensional configuration space 
with coordinates Xj. 

In Sec. 6, we present analogues of equations ( |2.2| ), ( |2.3| ) and ( p^.4| ) in the difference set- 
up. The isospectral flows are connected with eliptic solutions to the Volterra hierarchy. 



3 Double-Bloch eigenfunctions of the difference Lame 
operator and the spectral curve 

In this section we study Bloch eigenfunctions of the difference Lame operator. Following 
the general scheme outlined at the beginning of Sec. 2, we construct a family of eigenfunc- 



tions depending on E and two additional spectral parameters. All the three parameters 
are constrained by two equations which define the spectral curve. 

Consider the eigenvalue equation for the operator L ( |1.5| ): 

, OAx + iri)ei(x - (£ + l)r]) ,, , ^,, , , , 

The coefficient function is double-periodic. Therefore, it is natural to look for solutions in 
the class of double-Bloch functions [Q, i.e., such that iplx+l) = Biip{x), ip{x+T) = Bt-iIj{x) 
with some constants Bi, B^-. These are going to be the additional parameters Pi from the 
general scheme of Sec. 2. 

Consider the function 

Its monodromy properties in x are $(x + 1, C) = $(a;, (), $(x + r, (^) = e ^'^*^$(x, (), i.e., 
it is a double-Bloch function. Moreover, it is the simplest non-trivial (i.e., different from 
the exponential function) double-Bloch function since it has only one pole. This function 
serves as a building block for more general double-Bloch functions. 

Let £ be a positive integer. We employ the following double-Bloch ansatz for the i/j: 

^|J{x) = ir^/" Y. ^.(C, K, EMx - jr^, () , (3.3) 

i=i 

where (,K parametrize the Bloch multipliers of the function ijj{x): Bi = K~i , B^ = 
/^^g-27r«c Tiie coefficients Sj depend on the indicated parameters only. 

Substituting (|3.3|) into ( |3.1|) and computing the residues at the points x = jr], j = 
0, ...,£, we get i + 1 linear equations 

e 

Y,M,,s,=0, 1 = 0,1,. ..,£, (3.4) 

i=i 

for i unknowns Sj. Matrix elements Mj,- of this system are given by the formula 






M„ = K5,,i^^-E5^,,+K~^ .... . :'\^ . , k}H + 



+ ^ ^77^ a f \Q ff ■ — ■ I n N i(>i,o - di,i) . 3.5 



Here i = 0,1, . . . ,i, j = 1,1,..., L The overdetermined system ( p.4| ) has nontrivial 
solutions if and only if rank of the rectangular matrix Mjj is less than t. Let M'^^^ and 
M^^) be £ X £ matrices obtained from M by deleting the rows with z = and % = 1, 



respectively. Then the values of three parameters C, -^, -E for which eq. (|3.4|) has solutions 
of the form ( |3.3| ) are determined by the system of two equations: det M'^^^ = det M^^^ = 0. 
They indeed define a curve. 

To obtain the explicit form of these equations, we expand the determinants with respect 
to the first row. This yields an explicit characterization of the spectral curve summarized 
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in the theorem below. Hereafter, the "eUiptic factorial" and "elliptic binomial" notation 
is convenient: 

(3.6) 



n 
m 



[n]\ 



[m] \[n — m]\ 



Theorem 3.1 The difference Lame equation ^3. 1\ ) has douhle-Bloch solutions of the form 
(\3.3[ ) if and only if the spectral parameters (, K, E obey the equations 



j=0 



Af(E) = 0, 



(3.7) 



e+i 



J2{-iyK-^e^{c-3v)[3-i] 

j=0 



,W 



^7-iiiE)=0, 



iW 



where Aj (E) are polynomials of {£ — j)-th degree explicitly given by the determinant 
formula 

(hereO <s<i), A^ = 1. 



2i 
s 



det E5ij + YiTVx — ^^^d-^ 



\t+\-i\ 






ij+l 



(3.8) 



l<ij'<s 



Let us list some useful properties of the polynomials Al,-\E). First, they obey the 
recurrence relation: 



^^^i(^) = ^Mtm + #t41i(^) 



[2i-s] 



(3.9) 



iW, 



iW 



with the initial condition A),\E) = 1, A}Z^{E) = {[i]/[2i])E. Next, it is clear from j ^M 
that 

4li-E) = i~ir4tiE) , 0<s<i. (3.10) 



The equations ( |3.7| ) define a Riemann surface F, which covers the complex plane. The 
monodromy properties of the ^-function (see (^^) in the Appendix) make it clear that this 
surface is invariant under the transformation 

Ci — ^C + ^, K^ — >Ke^^'\ (3.11) 

The factor of the F over this transformation is an algebraic curve F, which is a ramified 



(|3.10|) that the curve 



covering of the elliptic curve with periods 1, r. It is clear from 
admits the involution 

{C,K,E)^^iC-K,-E), (3.12) 

so the spectrum is symmetric with respect to the reflection E — * —E. Another result 
of 1^, which is not so easy to see from ( |3.7| ), is that the curve F is at the same time a 
hyperelliptic curve. 



Theorem 3.2 The curve T is a hyperelliptic curve of genus g = 2i. The hyperelliptic 
involution is given by 

(C, K, E) ^ {2Nr^ ~C,K-\E), N = \i{i + 1) . (3.13) 

The points P = {(,K,E) & T of the curve parametrize double-Bloch solutions ^(x) = 
ijj{x, P) to eq. ( ^. 1\ ), and the solution iplx, P) corresponding to each point P G F zs unique 
up to a constant multiplier. 



For the proof see [Q. Here we give a few remarks. The involution ( 3.13|) looks best in 



terms of the function \l'(x) = ijj{x) Ilj=i 6i{x — jrj) which satisfies the eigenvalue equation 
L^ = E"^ with the L as in ( [1.3| ) (cf. ( |1.4| )). Then the hyperelliptic involution simply 
takes "^{x) to \E'(— a;). The genus of the curve can be found by counting the number of 
the fixed points of this involution. At the fixed points the two solutions, \l/(x) and \l/(— x), 
are linearly dependent: \&(— x) = r\l/(x). Writing out the eigenvalue equation at x = 0, 
we obtain the necessary condition \1'(— r/) = \l/(r/), so r = 1. The ansatz ( |3.3| ) for ip is 
equivalent to the ansatz 

^(x)=K^/^n^i(x+y,) 

i=i 
with J2j Vj = C- At K = 1, the dimension of the linear space of even functions of this form 
is known to be 2£ + 1. Adding images of the fixed points under the involution (|3.12|) , we 
eventually get 4£ + 2 fixed points, so, by the Riemann-Hurwitz formula, the genus is equal 
to 21 

Taking into account the symmetry E —>■ —E, we can write the equation of the hyper- 
elliptic curve in the standard form: 

2^+1 

w'= UiE'-E^). (3.14) 

The hyperelliptic involution takes {w,E) to {—w,E). 

In ( |3.14|) , w is the eigenvalue of a non-trivial operator W commuting with L on their 
common eigenf unction. The explicit form of the operator W was found in [pO|] : 






fc=0 



where ^i{x) = Uf=o Oi{x + (j - £)r/). 

Let us conclude this section by examining the behaviour of the spectral curve in the 
vicinity of its "infinite points", i.e., the points at which the function E has poles. From 
( p.7| ) we conclude that there are two such points: cxd+ = {( -^ 0, K -^ oo, E —^ cxd) and 
oo_ = (C ^ 2Nr], K ^0, E ^ oo). In the neighbourhood oi oo± E = K^^ + o{K^^). In 
terms of the variables (w, E) these points are oo± = (w ^- ±oo, E —>■ oo). 



4 Edges of bands 

The edges of bands ±i?j, i.e., the branch points of the two-sheet covering (|3.14|), are values 



of the function E = E{P) at the fixed points of the hyperelliptic involution. As is clear 
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from ( |3.13|) , the fixed points lie above the points ( = Nt] + Ua, where Ua are the half- 
periods: ui = 0, 1^2 = |, ^3 = |(1 + t), 1^4 = |t. The corresponding values of K are 
determined from ( |3.11| ). Then the set of the branch points Ei is fixed by Theorem 3.1. 



Corollary 4.1 Let Sa, a = 1, . . . ,A be the set of common roots of the polynomial equations 



j=0 



iW 



Af>{E)=0. 



(4.1) 



e+i 



Y.i-^yuiN-j)v)[j-i] 

j=0 



'+1 
J 



AtUE) = ^. 



where 6a are Jacobi 6-function^. Then the set of the edges of bands ±Ei is the union of 
Ua=i ^a and its image under the reflection E — > —E. 



Let us give two examples. At £ = 1, the set Si is empty while £a for a = 2, 3, 4 contains 
one point. From ( [4. ID we find: 



^ _o ^/3+l(^)^7+l(^) 

" ^^+i(0)^,+i(0) 



where {a, P, 7} is any cyclic permutation of {1, 2, 3}. At £ = 2, the set Si has two 
elements Ei, E2 obtained as solutions of the quadratic equation [2]E'^ + [2]^E + 2[4] = 0, 
so that 



''' 2 1^ 9U^) ^ efir^) %i{2r^) 
For a = 2, 3, 4 each set Sa has one element -Ea+i^ 

ei{2r])ea{27]) 



Ea+1 



Oliv)Oaiv) ■ 



In general, it is possible to prove 0] that 
i(£-l), £ odd 



#(^1 



i£+l, 



#(^2,3,4 



even 



Note that #(u^=i^a) = 2£ + 1 that agrees with ( pTTl) . 



Ui+l), ^ odd 



2 ^ ' 



even 



(4.2) 



5 Relation between the Bloch multipliers 



To simplify equations of the curve (|3.71 ), one can try to eliminate one of the variables and 
obtain a single equation for the other two. Here we show how to eliminate E. This leads 



^Definitions and transformation properties of the Jacobi ^-functions 9a{x\T), a = 1, 2, 3, 4, are listed in 
the Appendix. For brevity, we write 0a(x|r) = 9a(x). 



to a closed relation between the two Blocli multipliers of the function ( |3.3| ) (parametrized 
through ( and K). Its form (see (|5.6| ), (|5.^ ) below) suggests an interpretation in terms of 
hypothetical combinatorial identities for elliptic numbers. 

At the first glance, the elimination of E from eqs. ( |3.7|) is hardly possible at all. Never- 
theless, there is an alternative argument leading directly to the relation between the Bloch 
multipliers. Here it is more convenient to deal with the difference Lame operator in the 
gauge equivalent form ( |l.tj| ). Our construction is based on the following simple lemma [^. 

Lemma 5.1 Let \l/(x) be any solution to the equation 

in the class of entire functions on the complex plane of the variable x, then 

^{JV) = '^{-JV), J = 1,2,... J. (5.2) 

This assertion follows from the specific location of zeros and poles of the coefficients of 
eq. (|5.1|) . Indeed, putting a; = in ( ^.Ij) , we have ^{rj) = \l/(— r^). The proof can be 



completed by induction. At x = ±£77 one of the coefficients in the l.h.s. of (5J.) vanishes, 
so the chain of relations ( |5.2|) truncates at j = i. 

Remarkably, the conditions ( p.2| ) and the ansatz 

■^{x)=K^/'^il[e^{x-jr])] J2smiKXMx-mr],C) (5.3) 

\i=l / rn=l 



for \1/ (equivalent to the ansatz ( |3.3|) for ip) with the same function ^(x,z) given by ( p.2|) 
allow one to find the relation between the Bloch multipliers even without explicit use of 
the difference Lame equation (|5.1|) . Plugging ( |5.3| ) into (|5.2|), we obtain £ equalities (for 



m = l, 2,. ..,£): 

K-s^ = {-lYK-9,{2mv) ( {[ ^^T^^^^l E $(-(m + n)v, C) ^n • (5.4) 

This is a system of linear homogeneous equations for s„. It has nontrivial solutions if and 
only if its determinant is equal to zero, whence we obtain the equation connecting C, and 
K: 

det(ir2-(5„„ + G^„(C))^^ ^=0, (5.5) 

where 

G„n(C) = (-l)'+'[2m] ( n f^ ) -^(-(^ + ^)^> . 

This equation defines a curve, which is the image of the spectral curve F under the pro- 
jection that takes {(,K,E) to (C,-^)- 

The equation of the spectral curve ( p.5| ) can be represented in the form 

N 

E(-iycf (r^)0i(C - 2jr^)ir2(^-^) = 0, (5.6) 
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where N = ^i{i+ 1) and C- {i]) are some coefficients depending only on r] and r sucli tliat 

^] iv) = Cn-.j{^)-i Cq {t]) = 1. To see tliis, we expand tlie determinant ( |5.5| ) in powers 
of K witli tlie lielp of tlie identity 

/ e,{x, + X, + o \ er\0ei{c + 2T:^,Xi) - el{x,-x,) 

' I 9,{x, + xj) ^<^_^.<„ nr=i^i(2x,) Ihiix^ + xj) 

(a particular case of the formula for the elliptic Cauchy determinant). Let A be the set 
{1,2,..., i}. For any subset J C A, A\ J is its complement, and cr( J) = J2jej J- Evaluating 
the elliptic Cauchy determinant at Xn = —2rjn, we get: 

V / l<m,n<£ 

- e "-'^;':'""" /^— '(-!)"-' n n il^li. (5.7) 

JCA ^UU fceJfc'sAV N'^" '^ IJ 

Thus, the coefficient CJ reads 

cf- E n n ^. (5.8) 



The symmetry j *-^ N — j \s now transparent. Note that the sum in ( ^.8|) runs over 
partitions of the number j into distinct parts not exceeding L 

We note that for any elliptic module r 



.W 



iV 



.S^r= , (5.9) 



(the usual binomial coefficient). To see this, consider the limiting case r -^ zoo. We have: 
exp(— 7rir/4)6'i(a;|r) — ^ 2sin(7ra;) as r — > zcxo (see ([A2D), so 



sm(7rr7) 



Then (|5.9| ) follows from the identity [|21| 



JCA A;GJfc'GA\Jvl I''? l<j<k<e 

which is a specialization of the Crtype Weyl denominator formula^. A detailed combina- 
torial analysis of the limiting cases r -^ ioo or r ^ of the difference Lame operators can 
be found in [BlU. 



^I am grateful to A.N.Kirillov for pointing out this identity and drawing my attention to the paper 

ill. 
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6 Isospectral deformations of the difference Lame op- 
erator and locus equations 

Finally, let us comment on isospectral deformations of the difference Lame operator. We 
are going to present difference analogues of the operator ( |2.2| ) and of the locus equations 



Instead of ( p.l|) , consider the equation 

ip{x + T]) + c{x)ip{x — T]) = Eip{x) (6.1) 

with a more general coefficient c{x) , which is an elliptic function represented in the form 

c(x) - ^';;;;;g:;) "' . pw=T''.(x-x,). (6.2) 

Note that in the case of the difference Lame operator the configuration of zeros of the p{x) 
is very specific: 

p(x)= n 0,ix+ij + k-i-l)r]), (6.3) 

^<j<k<i 

SO all but two cancel in the c{x). 

The isospectral fiows are the fiows of the Volterra hierarchy for the c{x). The first 
equation of the hierarchy, 

OCiX} 

—^ = -c{x){c{x + r])~c{x-r])), (6.4) 

is the compatibility condition of the spectral problem ( |6.1| ) and the linear problem 

dipix) 



dt 



c{x)c{x — 'r])ip{x — 2?7) . 



Recall that changing the variables as t -^ \if't, x ^ x — 2rjt and setting c{x) = 1 — rfu{x), 
one gets the KdV equation for m as r^ — * 0. 



Substituting the pole ansatz (|6.2| ) into the Volterra equation and requiring the residues 
at the poles to be zero, we get the following two systems of equations (j = 1, 2, . . . , i£(£+l)): 



^(m)/2 



^;(o) . _'^'i^''e,{x,-Xk + 2T^)e,{xj-x,-r^) 



X 



n 



(6.5) 

e[{o) ^ _'(':^/^ e,{x,-Xk-2r])e,ix,-Xk + v) 

where Xj = dtXj. Solving these equations, or, equivalently, the Volterra equation with the 
initial condition ( |6.3D , one arrives at a family of isospectral deformations of the difference 
Lame operator. The two systems (|6.5|) must be satisfied simultaneously, therefore, the 
right hand sides are identical. Whence we obtain the necessary conditions for solutions to 
exist: 

'^'S^' e,{x,-xk + 2r^)el{x,-xu-r,) ^ ^ 

k}{^^j 9i{xj - Xfc - 27])ef{xj -Xk + v) 
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that is the difference analogues of ( p.3|) . The results of 0] imply that these equations 
define an equilibrium locus of the Ruijsenaars-Schneider model. Like in the differential 
case, the locus is not compact. Its closure contains, in particular, the point corresponding 
to the degenerate configuration ( |6.3|) . From the general arguments of ref. @] it follows 
that the connected component of the locus that contains this point at the boundary is 



^-dimensional. Expanding (|6.6| ) in r] ^ 0, we do get (2.3) 
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Appendix. Theta-functions 

We use the following definition of the Jacobi ^-functions: 

^M^) = -E exp i-^iTik + -f + 2m{x +-){k+ - 

^2(a;|r) = E exp iiiiTik + -f + li^ixik + -) j , 



1 



6*3 (x|r) = E exp [txit]^ -|- 2Tcixk] , 
fcez 

^4(a;|r) = ^ exp I nirk^ + 2Txi{x -\ — )k \ . 



keL 

They can be represented as infinite products 



oo 

2-ni(kT—x) 



^i(a;|r) = 2 sin(7rx)e"^"/^ 11 (^ - e^"'''") (^ - e2"*(^"+^)) (l - e 

k=\ 

oo 

^2(a;|r) = 2 cos(7ra;)e"^"/^n (l " e^"'''") (^ + e2"^('="+^)) (l + e^^^^^""^) 

k=\ 

oo 



fc=l 



0(J 



(Al) 



(A2) 
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The p-function is related to the 6i as follows: 

p{x\l/2, r/2) = -— log^i(x|r) + const . (A3) 

Throughout the paper we write 6'a(x|r) = 6'a(x). 

The transformation properties of the theta-functions for shifts by (half) periods are: 

Oaix ± 1) = (-l)^"^i+^"-20„(x) , Oaix ±t) = (-l)^'"-i+^°-4e-"^"^2™^<,(x) , (A4) 

ei{x ± f ) = ±ze-3'^*^^™^4(a;) , (A5) 
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